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Abstract 

Hyperfine splittings between the heavy vector {D*, B*) and pseudoscalar 
{D, B) mesons have been calculated numerically in lattice QCD, where the 
pion mass (which is related to the light quark mass) is much larger than its 
physical value. Naive linear chiral extrapolations of the lattice data to the 
physical mass of the pion lead to hyperfine splittings which are smaller than 
experimental data. In order to extrapolate these lattice data to the physical 
mass of the pion more reasonably, we apply the effective chiral perturbation 
theory for heavy mesons, which is invariant under chiral symmetry when the 
light quark masses go to zero and heavy quark symmetry when the heavy 
quark masses go to infinity. This leads to a phenomenological functional 
form with three parameters to extrapolate the lattice data. It is found that 
the extrapolated hyperfine splittings are even smaller than those obtained 
using linear extrapolation. We conclude that the source of the discrepancy 
between lattice data for hyperfine splittings and experiment must lie in non- 
chiral physics. 


PACS Numbers: 12.39.Fe, 12.39.Hg, 12.38.Gc, 12.40.Yx 


I. Introduction 


The past few years have seen much progress in lattice gauge theory, which is the 
only quantitative tool currently available to calculate nonperturbative phenomena in 
QCD from hrst principles. These phenomena include the spectrum of light hadrons 
[]^, the spectrum of heavy hadrons |^, the decay constants of heavy mesons 
[Q, nucleon structure functions and the Isgur-Wise function for B —> D{D*) 
transitions |P. 

In Ref. [^, the authors studied extensively the spectra of the D and B mesons 
using non-relativistic QCD (NRQCD) on the lattice in the quenched approxima¬ 
tion. For spin-independent splittings such as the splittings between strange and 
non-strange mesons, good agreement with experiments was obtained. However, it 
was found that the hyperhne splittings between D{B) and D*{B*) are much smaller 
than the experimental data. In fact, three lattice data for the mass of each of these 
mesons were obtained in the region where the mass of the pion is larger than about 
680MeV, which is much larger than the physical mass of the pion. With naive 
linear extrapolations from the unphysical region to the physical value of the pion 
mass, the extrapolated hyperhne splitting between D and D* is about llOMeV for 
(3 = 5.7 compared with the experimental value 140MeV, whereas for B and B* 
the extrapolated hyperhne splitting is about 29MeV compared with the experimen¬ 
tal value 46MeV. Obviously these large diherences between the extrapolated and 
experimental values merit careful investigation. 

It is well known that QCD possesses a chiral SU{3)l x SU{3)ii symmetry in the 
limit where the masses of light quarks m, d, and s go to zero. This symmetry is spon¬ 
taneously broken into SU{3)v, leading to eight Goldstone bosons. The interactions 
of these pseudoscalar mesons are described by an ehective chiral Lagrangian which 
is invariant under SU{3)l x SU{3)r. Another interesting quark mass limit in QCD 
is the heavy quark limit where the masses of heavy quarks, c and b, go to inhnity. 
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In this limit there are heavy quark flavor symmetry and heavy quark spin symmetry 
SU{2)f X SU{2)s. Based on these symmetries which are not manifest in the full the¬ 
ory of QCD, an effective theory for heavy quark interactions which is called heavy 
quark effective theory (HQET) was established 0. With the aid of HQET, the 
physical processes involving heavy quarks are greatly simplified. The interactions of 
heavy hadrons containing one heavy quark with the light pseudoscalar mesons vr, K 
and r] should be constrained by both chiral symmetry and heavy quark symmetry. 
The combination of these two symmetries leads to an effective chiral Lagrangian for 
heavy hadrons which is invariant under both SU (3)^ x SU (3)^ and SU (2) f x SU (2)^ 
transformations. There has been considerable work in this direction in recent years 


In the past few years there has been a series of work dealing with extrapolations of 
lattice data for hadron properties, such as mass |^, magnetic moments [^, parton 
distribution functions |l^, and charge radii |l^, to the physical region. In these 
extrapolations the inclusion of pion loops yields leading and next-to-leading non- 
analytic behavior. This leads to rapid variation at small pion masses while lattice 
data are extrapolated to the physical pion mass. However, when the pion mass 
is greater than some scale A, which characterizes the physical size of the hadrons 
which emit or absorb the pion, the hadron properties vary slowly and smoothly. It is 
obvious that extrapolations which ensure the correct non-analytic behavior of QCD 
in this way should be more reliable than a naive linear extrapolation. 

With these considerations in mind, the aim of the present work is to extrapolate 
the lattice data in Ref.[^ to the physical region, while building in the constraints 
of chiral perturbation theory for heavy mesons, and to compare the extrapolated 
hyperflne splittings with experiments. 

The remainder of this paper is organized as follows. In Section II we give a brief 
review for chiral perturbation theory for heavy mesons. In Section III we apply 
this theory to calculate self-energy contributions to heavy mesons from pion loops. 
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Based on this, we propose a phenomenological, functional form for extrapolating 
the lattice data for heavy meson masses to the physical region. Then in Section IV 
we use this form to £t the lattice data and give numerical results. Finally, Section 
V contains a summary and discussion. 

II. Chiral perturbation theory for heavy mesons 


When the heavy quark mass mg {Q = b, or c) is much larger than the QCD scale 
Aqcd, the light degrees of freedom in a heavy hadron are blind to the flavor and 
spin orientation of the heavy quark Q. Therefore, dynamics inside a heavy hadron 
remains unchanged under SU{2)f x SU{2)s transformations. In the opposite mass 
limit nig 0, the QCD Lagrangian possesses an SU{3)l x SU{3)r chiral symmetry. 
The light pseudo-Goldstone bosons associated with spontaneous breaking of chiral 
symmetry are incorporated in a 3 x 3 matrix 


S = exp 


/2fM' 


where is the pion decay constant, = 132MeV, and 


/ + 0, 


M = 


V2'' ' 

TT” 


TT 


K+ \ 




V 


K- 

Under SU{3)l x SU{3)ji transformations 


h) 


( 1 ) 


( 2 ) 




(3) 


where L G SU{3)l and R G SU{3)r. 

While discussing the interactions of Goldstone bosons with other matter fields 
it is convenient to introduce 

e = v^. (4) 

Under SU{3)l x SU{3)fi transformations 

e ^ = U^R^, (5) 
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where the unitary matrix U is a complicated nonlinear function of L, R, and the 
Goldstone helds. 

In order to discuss the interactions of Goldstone bosons with heavy mesons, 
which consist of a heavy quark Q and a light antiquark (a =1, 2, 3 for u, d, s 
quarks respectively), a 4 x 4 Dirac matrix Ha is introduced as follows: 

H.(v) = - ft75), (6) 

where P*^ and Pa are field operators which destroy vector and pseudoscalar heavy 
mesons with fixed four velocity v, respectively. P*^ satisfies the constraint 

v.Pr = 0. (7) 

Since Ha is composed of a heavy quark and a light antiquark, under SU (3 )l x 
SU{3)n 

Ha ^ ( 8 ) 

and under heavy quark spin symmetry 

Ha ^ SHa, (9) 

where S G SU{2)s. 

Dehning 

//a = 7oi/+7o, (10) 

we have 

Ha(v) = (p:;y + f’+75)^4, (11) 

which transforms as Ha —> UabHh and Ha HaS~^ under chiral symmetry and 
heavy quark symmetry, respectively. 

It is convenient to introduce a vector held 14^, 

= ( 12 ) 
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(13) 


and an axial-vector field 

A‘i, = \(ed‘‘i-id“eu. 

Under SU{?,)l x SU{V)r, V^^ f/U^[/+ + Ud>^U+, and UA>^U+. Defining 

the covariant derivative 

iD^H)a = d^Ha - (14) 

we find that {D^H)U+ nnder SU{3)l x SU{3)r. 

The Lagragian for the strong interactions of heavy mesons with Goldstone psen- 
doscalar bosons should be invariant under both chiral symmetry and heavy quark 
symmetry, since we are working in the limit where light quarks have zero mass and 
heavy quarks have infinite mass. It should also be invariant under Lorentz and par¬ 
ity transformations as required in general. The most general form for the Lagragian 
satisfying these requirements is [§] 

C = -TiiHaiv^{D^H)a] + gTT{HaH,^^Al-f,), (15) 


where g is the coupling constant describing the interactions between heavy mesons 
and Goldstone bosons. Obviously, g is universal for D, B, D*, and B*. Since g 
contains information about the interactions at the quark and gluon level, it cannot 
be fixed from chiral perturbation theory for heavy mesons, but should be determined 
by experiments. 


From Eq.(|T^, the propagator for the pseudoscalar meson D (or B) is 


2v ■ p 


where p is the residual momentum of the meson. For the vector meson D* (or B*), 
the propagator is 


- Vf^Vu) 

2v ■ p 

In the limit mg —> oo, there is no mass difference between pseudoscalar and vector 
mesons. 
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Since in this work we will study hyperfine splittings, we need to take l/rriQ 
corrections into account. At order l/mq in HQET, the term which is responsible 
for hyperfine splittings is the color-magnetic-moment operator, (where 

hy is the heavy quark field operator in HQET and is the gluon field strength 
tensor). This leads to the following correction term to C in Eq.(^): 

(16) 


where A 2 is a constant which also contains interaction information at the quark and 
gluon level, and which is same for D, B, D*, and B* at the tree level. (When QCD 
loop corrections are included, A 2 depends on mg logarithmically.) Finally, we note 
that the inclusion of the other term at order l/mg in HQET, leads 

to a slight mg dependence of the coupling constant g. 

Adding the term ([IBD to the Lagrangian for HQET, Eq.([l5|), and using Eqs.(|) 
and ([Tl|) , we can easily see that the mass difference between vector and pseudoscalar 
heavy mesons is 

A 

mg 

and consequently, the propagators for heavy pseudoscalar and vector mesons become 

i 

2(n - pT |A) 



-ijOfiu - Vi^Vu) 
2{vp-\A) ’ 


respectively. 

In order to consider the interactions of heavy mesons with Goldstone bosons, we 
substitute ^ = exp{iM/f.^) into Eqs.(|^, |^) and obtain the following expressions 
for and A^: 


V, = ^[M,d,M]+0{M^), 

TV 


( 18 ) 

( 19 ) 
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Substituting Eq.([T8|) and Eq.([T^) into Eq.(p!5D we have the following explicit form 
for the interactions of heavy mesons with Goldstone bosons: 

S„M|j.(P;+Pr - P+Pj) 

J TT 

In 

+i£»“''”pyPj;«„a„Afj„), (20) 

where 0{M^) terms are ignored. 

Chiral symmetry can also be broken explicitly by nonzero light quark masses. 
Under SU{3)l x SU{3)ji, light quark mass terms transform as (3^, 3k) + (3k, 3k). 
Then, to leading order in the explicit chiral symmetry breaking from light quark 
masses, the following terms are added to the Lagrangian in Eq. (p!5D : 

XiTiHbHai^mq^ + + X^TiHaHai^mq^ + ( 21 ) 

where Ai and A'^ are parameters which are also independent of the heavy quark mass 
in the limit mg —>■ 00 . 

III. Formulas for the extrapolation of heavy meson masses 

In this section we use the Lagrangian for the interactions of heavy mesons with 
light Goldstone bosons to calculate pion loop corrections to the masses of heavy 
vector and pseudoscalar mesons. This yields the dependence of heavy meson masses 
on the pion mass. Then we propose a phenomenological functional form for extrap¬ 
olating lattice data for heavy meson masses to the physical pion mass. 

From Eq.(^) we can see that there are hve possible diagrams for pion loop 
corrections to heavy meson masses. These diagrams are shown in Fig.l for D and 
B mesons, and in Fig.2 for D* and B* mesons. Fig.1(a) and Fig.2(a) arise from the 
first term in Eq. (^0|) . It can be easily seen that these two diagrams do not contribute 
to the masses of heavy mesons and we will not consider them from now on. 
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Fig. 1(b) represents the pion loop correction to the heavy pseudoscalar meson 
propagator (P could be P or P). In momentum space it can be expressed as 


2(!i-p+ |A) 


;-2iE) 


2(ti ■p + |A)’ 


( 22 ) 


where p is the residual momentum of the pseudoscalar meson and —is given by 
the following integral: 

3^2 r d^k k^-{vkf 


- 2iE = 


(23) 


fi J (2!r)‘ Iv-(p-k)- jA](F - mj) ’ 
where k is the momentum of the pion in the loop, and is the pion mass which 
is not necessarily its physical mass. It can be seen that S is a function of n • p, 
and A. After the correction from S is added, the heavy meson propagator is 
proportional to 

--T- P'l) 

V ■ p — mo — ^[v ■ p) 

where mo is the mass term without S correction (for the propagator of the pseu¬ 
doscalar heavy mesons, mo = — |A). The physical mass of the heavy meson, m, is 
defined by 

[vp-mo-T.{v ■p)]\y.p=m = 0- (25) 

Therefore, to order 0{g‘^) we have 

m = mo -l- S(n ■ p = mo). (26) 

In order to calculate the integral in Eq.(^3l), we need to deal with the following 
integral: 


d^k 


k^k’' 


~ J (27r)4 [vk-6] {k^ - ml) ’ 

where 6 = v ■ p — |A for Eq. (^3|) . On the grounds of Lorentz invariance, in general 
we have 

+ X2V^^v^, (28) 


where Xi and X 2 are Lorentz scalars, which are functions oiv-p, ml, and A. It can 
be easily seen from Eq.(|23|) that S ~ {g^'^ — v^v’^)Xp,^. Consequently, only the Xi 








term contributes, since = 1. Multiplying both sides of Eq.(^) with — v^v 
we have 

k'^ — {v ■ k)‘^ 


^■‘ = 3 


d^k 

{2Tr)*v^{ko 


v-kH-(5 

7,0 


ie){k-^ - Wi + ie)' 


(29) 


where W? = IkP + 


mz 


Next we carry out the integration over fcg by choosing the appropriate contour. 
From Eq.(^), when we calculate mass corrections from pion loops, we need the 
value of S at n • p = niQ. Since n ■ p is a Lorentz scalar, we are free to choose some 
special value of v for this purpose. With = 1 and v = 0, and choosing the contour 
in the lower half plane, in which there is only one pole for feg, Wk — we arrive at 
the following expression: 

'>• r d^k IkP 


= 


Qj (27r)3Wfc(Wfc-<5)- 


(30) 


In Refs. []^|-1^1, it has been argued that when the Compton wavelength of the 
pion is smaller than the source of the pion, pion loop contributions are suppressed as 
powers of m.^/A where A characterizes the hnite size of the source of the pion. We 
follow this argument and introduce a cutoff A in the integral Since the leading 
non-analytic contribution of pion loops is associated with the infrared behavior of 
the integral, it does not depend on the details of the cutoff. In the following, we will 
treat A as a parameter to be hxed by lattice data. 

With the cutoff A, we obtain the following result for the integral (^Ol): 

— 6 


Xi = 


12{ml - 


727r2 ^ 

+3(5(25^ — 3m^)ln 


A + JA'^ + ml — 6 

arctg- *, =-arctg 


A + wA^ + 




^Jml - ( 5 ^ ^Jml - 6^ 

+ 3(5A\/A2 + ml + 6(3=“ - ml)A + 2Ad , 


(31) 


when ml > S^; 


m. 


()3/2i 


n 


'A ++ ml-5-^ 

A + ^A^ + ml-5+^ 


mt 


mt 




6 + \J^‘^ — ml 


m^ 


5 — \J^‘^ — ml 
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\ _|_ /A2 _|_ 777,2 _ 'I 

+3(5(25^ — 3m^)ln- - - - + 3SAJA'^ + ml + 6(5^ — ml)A + 2A^ > , 


when 

In the case where 5 = 0, we have 


^1 = 


(3m^arctg— - 3m^A + A' 
3o7r"' V 


(33) 


If we take the chiral limit m„ 
non-analytic term is 


0 in Eq.(|3^), we can see that the leading 


(34) 


? m"* 


327r2 5 

It can be easily checked that the same chiral limit behavior is obtained if we work 
with the dimensional regularization method in evaluating This is because the 
leading non-analytic contribution of the pion loops is only associated with their 
infrared behavior. 

In Eg. (^Bf) , mo for a pseudoscalar heavy meson is — |A, and hence 6 in Eq.(^) 
equals —A. Therfore, for pion loop contributions to the mass of a pseudoscalar 
heavy meson P, we have 

m^ A 


CTp = 


g 

167r2/2 


12{ml - A2)3/2 


A -h aA2 + ml + a 

arctg-=-arctg 

m2 - A2 


ml — A2 


-3A(2A2 - 3m2)ln 

+2A^} , 


A + wA2 ml 


m. 


3AAy^A2 -1- ml -I- 6(A^ — m^)A 


(35) 


when m^ > A^; 


CTp = 


In 


'A + ^A^ + ml + A- ^ 


mt 


m. 


+ A + JA2 


mt 


m^ A - 


167r2/2 \^A + -I- m2 -f- A -I- y^A2 — m' 

yY -|- /y^2 I ‘fffp' 

x6(A^ — — 3A(2A^ — 3m^)ln- - - - — 3AAXa2 -|- ml 

m^ ^ 

+Q{A^ -ml)A+ 2A^], (36) 


when ml < A^. 


( 32 ) 
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Now we turn to the pion loop corrections to heavy vector meson masses. First 
we discuss Fig.2(b), where P could be D or B. This diagram is caused by the PP*7i 
vertices in Eq.@). In momentum space, Fig.2(b) can be expressed as 


-KOlip - Vf^Vp) -i{gau - VaVu) 


-{2iW^)- 


(37) 


2{vp-\A) 2{vp-\A) ’ 

where p is the residual momentum of the heavy vector meson and 2in^”' is given by 
the following integral: 


2 m- 


(38) 


/2 J (27r)4 [v.{p-k) + I A](fc2 - ml) ’ 
with k the momentum of the pion in the loop. While we evaluate Eg. (|3^) , 5 in 


Eq.(^) is equal to v ■ p + |A. Because of the factor Qpp 


Vf_,Vp 


in Eq. ( 0 ). only the 


Xi term contributes. We dehne the coefficient of in 11^'^ to be fl 


= n, 


(39) 


SO that after Fig.2(b) is included, the propagator of a heavy vector meson becomes 
proportional to Eq. (P^) , with S being replaced by If and mo being equal to |A. 
While calculating the pion loop corrections to a heavy vector meson, we have to £x 
V - p = I A, as required in Eq.(^). Consequently, 6 in Eq. (^7D is equal to A. Then 
with the aid of Eqs.(^ |3^ we have 

' A 


n = 




I2{ml - A‘^)^A 


arctg- 


A + + ml 


mt 


A2 


arctg 


m^ 


A 


mt 


A2 


yY -|- I _ 'I 

+3A(2A2 - 3m2)ln-^ + 3AA Ja^ + ml + 6(A2 - ml)A + 2AH , 

m^ '' J 

(40) 


when ml > A^; 

,2 


n = - 


g 




In 


'A + ^A^ 


mt 


A 


A2 


mt 


A + AA^ + m2 — A + JA^ — m: 


m^ - A + AA2 


mt 


m-n- — A — s! X? — m2 


x6(A^ — ml)^A _|_ 3A(2A^ — 3m^)ln- ^ —- + 3AAi/A2 + ml 

m^ ^ 

+6(A2-m2)A + 2A=^}, 


(41) 
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when < A^. 

Now we discuss Fig.2(c), which arises from the P*P*7r vertices in the Lagrangian 
(1^) . In momentum space, Fig.2(c) can be expressed in the following explicit form: 

-Vf,Vp,) f3g^ 


2(vp- iA) \ /2 


aiZ/ipicri Q2^'2P20-2/^ _ \ 

^ ^ \y<Jip2 ^ai ^p2 J ^1^1 


d^k 


h h 
^0.2 


(42) 


(43) 


'^i9l'U2 '^vV52) 

J {27i)^ [v ■ {p - k) - \A]{k^ - ml) j 2{v ■ p - \A) 

After some algebra, this expression can be written in the form: 

-i{9pu - VpvP) i ^ 

2(n-p-iA) 2(n-p-iA) /2 

where m. Xi, 5 = v ■ p — |A. 

Because of Fig.2(c), the propagator for a heavy vector meson becomes propor¬ 
tional to 

TT fj7’ 

V ■ p — IA — If' 

where 

(45) 


n' = 

J TT 


and again 5 = v ■ p — |A in Xi. When we calculate corrections to a heavy vector 
meson mass from Fig.2(c), v ■ p should be taken to be |A, as required by Eq.( 
Hence 5 = 0. Using Eq. (|^ we have 

A 


n' = 


9 

'I2vr2/^ 


SmtaTctg- - 3m k + A' 

m^ 


(46) 


Adding H in Eqs. (^0]) and (|4l|) and H' in Eq. (^6[) together, we have the following 
expression for pion loop contributions to the mass of a heavy vector meson P*: 


dp* — 


167r2/2 


4{ml - A2)3/2 


A + jA‘^ + ml-A m^-A 

arctg- , =-arctg- 


mt — A^ 


ml — A^ 


A r, n 2\ -|- 4 /A 2 -|- ml , - 

-|-4m^arctg-h A(2A^ — 3m^)ln- - -h AAJA^ -|- ml 

m^ m^ * 


+2(A2 - 3ml)A + 2A3} , 


(47) 
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when > A^; 


(T p* — 




In 


A + a/A^ + rni — A — -i/A^ — m- 


A + a/A^ + 772,2 — A + a/A2 — 772 


772jr — A + a/A2 — 772 


772,^ — A — a/A2 — 772 


A A. -1- 1 -l- T?7^ 

x2(A^ — 772^)^^^ + 4772^arctg-h A(2A^ — 3772^)ln- - -^ 

727 ^ 727 ^ 

+AA^A2 + m2 + 2(A2 - 3m2)A + 2A3| , (48) 

when 727^ < A^. 

Eqs.(^ ^ ^ are valid provided 772 ,^ is not far away from the chiral limit, 
i.e., when 772 ,^ < A. It can be seen from Eq.(p9D that pion loop contributions vanish 


in the limit 772 ^ —> cx). In order to extrapolate lattice data from large pion mass to 
the physical mass of the pion we need to know the behavior of heavy meson masses 
at large 772 ^. As we know, a heavy meson is composed of a heavy quark and a light 
quark. In HQET, the heavy quark mass, mg, is removed. Therefore, as the light 
quark mass, rriq, becomes large, the heavy meson mass should be proportional to mg. 
On the other hand, explicit lattice calculations show that over the range of masses of 
interest to us, 772 ^ is proportional to rriq |^. Consequently, as m.^ becomes large, the 
heavy meson mass becomes proportional to 772 ^. Based on these considerations, in 
order to extrapolate lattice data at large ttTtt to the physical value of the pion mass, 
we propose the following phenomenological functional form for the dependence of 
heavy meson masses on the mass of the pion over the mass range of interest to us: 


772p = ap + bpml + ap, 


(49) 


for heavy pseudoscalar mesons, and 

mp* = ap* + bp*ml. + ap*, (50) 

for heavy vector mesons, ap and ap* are given in Eqs.(|3^, |36|) and Eqs.(|47|, |^), 
respectively. In next section, we will use this form to extrapolate lattice data for 
heavy mesons. 
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Before turning to the lattice data we comment briefly on the explicit chiral 
symmetry breaking by the terms in Eg. (|2l|) . Substituting Eqs.(|^, |^) into 

Eq.(^) we have the following explicit expression for Eq.( 


4A. f: - p:Pa) +4a; •£ hp^p:" - ( 51 ) 

a=l a=l a=l 

where we have made a Taylor expansion for ^ and omitted 0(1//^) terms. It can 
be seen clearly that Eq. contributes equally to the mass of a heavy pseudoscalar 
meson and that of a heavy vector meson. Therefore, it does not contribute to their 
mass difference. Corrections to Eq.(|^) may modify the propagators of heavy mesons 
to order mgO{l/f^), with extra suppression from mg with respect to the pion loop 


effects in Eqs. 


and Eqs.(^, and we will ignore them. 


IV. Extrapolation of lattice data for heavy meson masses 


The spectra of the D and B mesons were calculated in Ref. [§], where NRQCD 
was used to treat the heavy quarks. In fact, when the mass of a heavy quark, mg, 
is much larger than Aqcd, h becomes an irrelevant scale for the dynamics inside a 
heavy hadron. Consequently, heavy meson states can be simulated on lattices with 
the aid of NRQCD (where mq is removed from the Hamiltonian) when the lattice 
spacing is larger than the Compton wavelength of the heavy quark. Two values 
for the bare gauge coupling /3, 5.7 and 6.2, were used. We choose to use the data 
at = 5.7, since in this case the inverse lattice spacing a~^ is about 1.116GeV, 
which is smaller than the masses of either the 6 or c quarks. The box size is 2.1fm, 
corresponding to the volume 12^ x 24. In their simulations, three different values for 
the hopping parameter k, 0.1380, 0.1390, and 0.1400, were used. The light quark 
mass is related to k as = ^(I/k — 1/Kc), with Kc = 0.1434. 

In our model there are three parameters to be hxed, ap(p*), bp(^p*), and A in 
Eqs. (^91) and (|^). These parameters are related to g and A 2 which represent in¬ 
teractions at the quark and gluon level and cannot be determined from chiral per- 
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turbation theory for heavy mesons. As discussed in Section II, they may depend 
slightly on mq. In our htting procedure, we treat them as effective parameters and 
assume that their mq dependence has been considered effectively. From the upper 
limit for the experimental decay width of D* —>■ Dn [Q we have the upper limit 
< 0.5 P, |T^. In our numerical work we let vary from 0.3 to 0.5. From Eq.(|T7|) 
A 2 is related to the mass splitting between a heavy vector meson and a heavy pseu¬ 
doscalar meson. From experimental data for the case of B mesons, where l/rriq 
corrections can be safely neglected, the value of A 2 should be around -O.OSGeV^. To 
see the dependence on A 2 of our hts, we allow A 2 to vary between -O.OSGeV^ and 
-0.02GeV2. 

Using the three simulation masses for D, D*, B, and B* in Ref.[^, we £x the 
three parameters, ap(p*), bp(^p*), and A in Eqs.(|i9|) and (|50|) with the least squares 
fitting method. In Tables 1 and 2, for different values of A 2 and g"^, we list the 
parameters ap(^p*), &p(p*), and A obtained by htting the lattice data. Furthermore, 
we give the masses of D, D*, B, and B*, and hyperhne splittings at the physical 
pion mass, after extrapolation. With these parameters ap(^p*), hp(^p*), and 

A, we obtain the masses of D, D*, B, and B* as a function of the pion mass, for 
different values of A 2 and g. These results are shown in Figs.3 and 4. The mass 
differences between D (B) and D* {B*) as a function of the pion mass are shown in 
Figs.5 and 6. From these hgures we can see that when the pion mass is smaller than 
about 600MeV the extrapolations deviate signihcantly from linear behavior. This 
is because the pion loop corrections begin to affect the extrapolations around this 
point. As the pion mass becomes smaller and smaller pion loop corrections become 
more and more important. Gonsequently, the dependence of these plots on A 2 and 
g^ is stronger when the pion mass is smaller. 

As discussed in Section III, the parameter A characterizes the size of the source 
of the pion. Since the sizes of D, D*, B, and B* are different, we expect A for them 
are different. This is consistent with what is shown in Tables 1 and 2. Furthermore, 
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we can see that the A difference between D and D* is mnch bigger than that between 
B and B*. This is because the size difference between a heavy pseudoscalar meson 
and a heavy vector meson is caused by l/mg effects. 

In the naive linear extrapolations pion loop corrections are ignored. Hence the 
results do not depend on A 2 and g^. In Table 3 we list the results of the linear 
extrapolations for comparison. 

It can be seen clearly that the mass differences between D [B) and D* (B*) in our 
model are even smaller than those obtained in the naive linear extrapolations. Since 
in the linear extrapolations, the hyperfine splittings at the physical mass of the pion 
for D and B mesons are already smaller than experimental data, the inclusion of 
pion loop effects makes the situation even worse. As shown in Table 1, the hyperhne 
splitting between D and D* is 0.080 ~ O.OOlGeV in the range of our parameters, 
compared with the experimental value 0.14GeV, while the result from the linear 
extrapolation is 0.114GeV in Table 3. Similarly, the hyperfine splitting between B 
and B* in Table 2 is 0.023 ~ 0.026GeV in the range of our parameters, compared 
with the experimental value 0.046GeV and the result from the linear extrapolation, 
0.031GeV in Table 3. 

In our fits, in addition to the uncertainties in the three parameters, ap(p.), 6p(p*), 
and A, which are caused by the uncertainties of A 2 and g, the errors in the lattice 
data can also lead to some uncertainties in these three parameters. Since all the 
three lattice data are at large pion masses, and since A is mainly related to the data 
at small pion masses, the error in our determination of A is very large. Besides A, 
ap(p*) and 6p(p*) also have some errors. In fact, the numerical results for ap(p*), 
6p(P*), and A in Tables 1, 2, and 3 correspond to the central values of the lattice 
data. As a consequence of heavy quark symmetry, the dynamics inside a heavy 
meson does not depend on the mass of the heavy quark, so we expect that the 
values of A for D and B mesons are not very different from those of light mesons 


T0| . Gonsequently, we expect that the values of A in Tables 1 and 2 are quite 
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Table 1: Fitted parameters, masses of D and D* and their difference at 
(Numbers without (with) brackets are for D {D*)) 



A2=-0.02GeV2 

A2=-0.02GeV2 

A2=-0.03GeV2 

A2=-0.03GeV2 


= 0.3 

^2 = 0.5 

= 0.3 

c/2 = 0.5 

A(GeV) 

0.43 [0.65] 

0.38 [0.55] 

0.45 [0.63] 

0.39 [0.55] 

flD[D*](GeV) 

0.5439 [0.6848] 

0.5436 [0.6796] 

0.5444 [0.6821] 

0.5438 [0.6807] 

bn[D*]{GeV ^) 

0.2082 [0.1783] 

0.2083 [0.1809] 

0.2078 [0.1798] 

0.2082 [0.1800] 

rn%D*] (GeV) 

0.5374 [0.6282] 

0.5364 [0.6233] 

0.5375 [0.6254] 

0.5369 [0.6166] 

mp. — rrip (GeV) 

0.0908 

0.0869 

0.0879 

0.0796 


Table 2: Fitted parameters, masses of B and B* and their difference at (Num¬ 

bers without (with) brackets are for B {B*)) 




A2=-0.02GeV2 

A2=-0.02GeV2 

A2=-0.03GeV2 

A2=-0.03GeV2 



g^ = 0.3 

g^ = 0.5 

g^ = 0.3 

g^ = 0.5 

A(Ge 

V) 

0.62 [0.65] 

0.53 [0.56] 

0.62 [0.65] 

0.54 [0.56] 

ClB[B* 

l(GeV) 

0.7540 [0.7915] 

0.7507 [0.7884] 

0.7535 [0.7918] 

0.7519 [0.7887] 

bB[B*] 

(GeV^) 

0.1657 [0.1581] 

0.1673 [0.1593] 

0.1661 [0.1579] 

0.1665 [0.1591] 

m%B 

*l(GeV) 

0.7138 [0.7394] 

0.7107 [0.7342] 

0.7146 [0.7389] 

0.7107 [0.7336] 

mp* • 

-m^*(GeV) 

0.0256 

0.0236 

0.0242 

0.0229 


reasonable. In our fits, we find that the errors of the lattice data cause 1% ~ 2% 
relative uncertainties for ap(p*) {P = D or B), 6% for bn, 9% for bn*, 8% for bn, 
and 9% for bn*- These lead to about 1.3% uncertainties for mn, 1.5% for m*n, 1.2% 
for mn, and 1.2% for mp*. As a result, the relative uncertainties of the hyperfine 
splittings are about 13% and 54% in the cases of D and B mesons, respectively. In 
spite of all these errors, our conclusion that the hyperfine splitting obtained after a 
careful treatment of chiral corrections are smaller than those obtained using naive 
linear extrapolation, remains unchanged. 
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Table 3: Fitted parameters, masses of P and P* and their difference at for 

linear extrapolation (Numbers without (with) brackets are for P {P*)) 


aD[D*' 

l(GeV) 

bDiD*]{GeV-^) 


m% —171^ {GeV) 

0.5397 

[0.6540] 

0.2112 [0.1995] 

0.5438 [0.6579] 

0.1141 


(GeV) 

&B[B*](GeV 


m% — m^{GeV) 

0.7311 

[0.7621] 

0.1812 [0.1780] 

0.7346 [0.7656] 

0.0310 


V. Summary and discussion 

QCD possesses chiral symmetry when light quark masses go to zero and heavy 
quark symmetry when heavy quark masses go to inhnity. Combining these two 
symmetries leads to chiral perturbation theory for heavy mesons which are invariant 
under both chiral symmetry and heavy quark symmetry. We have evaluated pion 
loop corrections to heavy meson propagators with the aid of chiral perturbation 
theory for heavy mesons as the Compton wavelength of the pion becomes larger than 
the size of the heavy mesons. This leads to the dependence of heavy meson masses on 
the pion mass. In order to study hyperhne splittings, we took the color-magnetic- 
moment operator at order l/mq in HQET into account. This operator breaks 
heavy quark spin symmetry and is primarily responsible for the mass difference 
between a heavy pseudoscalar meson and a heavy vector meson. The small masses 
of the light quarks break chiral symmetry explicitly. We showed that contributions 
to the mass difference between a heavy pseudoscalar meson and a heavy vector 
meson from these terms are suppressed by light quark masses with respect to the 
pion loop contributions we considered in chiral perturbation theory. When 
becomes large, lattice data show that heavy meson masses are proportional to 
Based on these considerations, we proposed a phenomenological functional form 
with three parameters to extrapolate the lattice data. Because it guarantees the 
model independent chiral behavior of QCD, our model is more appropriate than 


18 



a naive linear extrapolation. The parameters in our model are fixed by the least 
squares fitting method, while fitting the lattice data for the masses of heavy mesons 
D, D*, B, and B* in the large pion mass region (m.^- > 680MeV). It is found that the 
hyperfine splittings extrapolated in this way are even smaller than those obtained 
in the linear extrapolations, in which the extrapolated hyperfine splittings for both 
D and B mesons are already smaller than the experimental data. 

There are some uncertainties in our model. We have two parameters, A 2 and g, 
which are related to the color-magnetic-moment operator at order l/mq in HQET 
and the interactions between heavy mesons and Goldstone bosons in chiral perturba¬ 
tion theory, respectively. In the ranges we choose for these two parameters, we have 
about 13% and 11% uncertainties for the hyperfine splitting in the D and B cases, 
respectively. Furthermore, the errors in the lattice data also lead to some uncer¬ 
tainties when we fix the three parameters ap(p*), &p(p*), and A in our model. Since 
all the lattice data are at high pion masses, the error in A is very large. However, 
we believe that the range of A we obtained is appropriate because of considerations 
based on heavy quark symmetry. The errors for ap(p*) and &p(p*) lead to about 
13% and 54% uncertainties for hyperfine splittings in the cases of D and B mesons, 
respectively. Despite all these uncertainties, the hyperfine splittings obtained in 
our model are smaller than those in the naive linear extrapolations. Our analysis 
shows that the current lattice data for hyperfine splittings at large pion masses are 
probably too small to give hyperfine splittings at the physical pion mass which are 
consistent with experiments. 

Some approximations made in current lattice simulations may be the cause of 
these small hyperfine splittings. The quenched approximation might be one reason. 
In fact, in lattice simulations for light spectroscopy the hyperfine splittings are 
also too small [^. Furthermore, as pointed out in Ref. p, the lattice results for 
hyperfine splittings are sensitive to the coefficient of the a ■ B term in NRQCD which 
is at order l/mq and which is the leading term to cause hyperfine splittings. The 
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inclusion of radiative corrections beyond tadpole improvement for this coefficient 
may increase hyperfine splittings. Another reason might be the light quark mass 
dependence of the clover coefficient in the clover action for light quarks. In addition, 
hnite size effects and higher order terms in NRQCD in lattice simulations may lead 
to an underestimate of hyperhne splittings at large light quark mass as well. More 
careful lattice simulations with more data and better accuracy are urgently needed 
to resolve this important problem. 
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Figure Captions 


Fig.l Pion loop corrections to the propagator of heavy pseudoscalar meson P {P 
could be D or 5). 

Fig.2 Pion loop corrections to the propagator of heavy vector meson P* [P could 
be D or B). 

Fig.3 Phenomenological hts to the lattice data for the masses of D* (the upper 
lines) and D (the lower lines) as a function of the pion mass. The solid (dashed) 
lines correspond to A 2 = —0.02GeV^ and = 0.3 = 0.5). The dot (dot dashed) 

lines correspond to A 2 = —0.03GeV^ and g^ = 0.3 {g"^ = 0.5). 

Fig.4 Phenomenological hts to the lattice data for the masses of B* (the upper 
lines) and B (the lower lines) as a function of the pion mass. The solid (dashed) 
lines correspond to A 2 = —0.02GeV^ and g^ = 0.3 {g"^ = 0.5). The dot (dot dashed) 

lines correspond to A 2 = —0.03GeV^ and g^ = 0.3 {g"^ = 0.5). 

Fig.5 Phenomenological hts to the lattice data for the hyperhne splitting between 
D* and iA as a function of the pion mass. The solid (dashed) lines correspond to 
A 2 = —0.02GeV^ and g^ = 0.3 {g"^ = 0.5). The dot (dot dashed) lines correspond to 

A 2 = —0.03GeV^ and g"^ = 0.3 {g"^ = 0.5). 

Fig.6 Phenomenological hts to the lattice data for the hyperhne splitting between 
B* and i? as a function of the pion mass. The solid (dashed) lines correspond to 
A 2 = —0.02GeV^ and g^ = 0.3 {g"^ = 0.5). The dot (dot dashed) lines correspond to 

A 2 = —0.03GeV^ and g"^ = 0.3 {g'^ = 0.5). 
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